Nonlinear oscillator with parametric colored noise: some analytical results 
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The asymptotic behavior of a nonlinear oscillator subject to a multiplicative Ornstein-Uhlenbeck 
noise is investigated. When the dynamics is expressed in terms of energy-angle coordinates, it is 
observed that the angle is a fast variable as compared to the energy. Thus, an effective stochastic 
■ dynamics for the energy can be derived if the angular variable is averaged out. However, the 

standard elimination procedure, performed earlier for a Gaussian white noise, fails when the noise is 
colored because of correlations between the noise and the fast angular variable. We develop here a 
specific averaging scheme that retains these correlations. This allows us to calculate the probability 
distribution function (P.D.F.) of the system and to derive the behavior of physical observables in 
the long time limit. 
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I. INTRODUCTION 



The mechanical model of a particle trapped in a nonlinear confining potential and subject to external noise has been 
widely studied to illustrate the interplay between randomness and nonlinearity in a classical dynamical system 
Because of the external noise, some parameters of the system fluctuate with time and the equation of motion becomes 
a stochastic differential equation with multiplicative (i.e., parametric) noise. Models with multiplicative noise can 
undergo purely noise- induced phase transitions 0,0,0 that can be realized experimentally, for example, in nonlinear 
electronic circuits Q or in systems with hydrodynamic instabilities : 8|. Besides, the presence of multiplicative noise 
provides a microscopic model for the anomalous diffusion of a particle in a fluctuating field of force in the limit of 
vanishing damping rate 0, ^3 ■ 

When the external randomness is represented by a Gaussian white noise, many analytical results can be derived 
from the Fokkcr-Planck equation. For example, in [TT| . we have calculated the growth exponents and the associated 
generalized diffusion constants of a nonlinear oscillator with random frequency in the limit of vanishing damping 
rate. The key feature of the method is to derive an effective first order Langevin equation for the action variable 
by averaging out the fast angular motion. The effective low dimensional problem can then be solved analytically. 
However, this method fails for colored noise. In fact, when the noise is colored, a new time scale, the coherence time 
, of the noise, appears in the problem and straightforward averaging over the fast angular variable leads to erroneous 
results because correlation terms between the noise and the fast variable are eliminated. We could analyze only 
qualitatively ^2] through self-consistent scaling arguments the long time behavior of a nonlinear oscillator subject to 
colored parametric noise and found that it is radically different from that observed with white noise. 

In this work, we develop a consistent adiabatic averaging procedure that allows us to derive analytical results for the 
general nonlinear oscillator subject to an Ornstein-Uhlenbeck parametric noise. We calculate, in the long time limit, 
the Probability Distribution Function (P.D.F.) and derive explicit formulae for the mean values of physical observables 
such as the energy, the velocity-square, the amplitude-square. The method we describe here is a generalization of the 
averaging technique introduced in [lflj for the classical pendulum with fluctuating frequency. 

The outline of this paper is as follows. In section 2, we describe the model under study and review results that 
O ■ have been obtained in previous works : we write the dynamical equation in terms of energy and angle coordinates and 
observe that a separation between fast and slow variables occurs in the long time limit. We then give the analytical 
• »h , expression of the P.D.F. and the scaling of physical observables in the white noise case. In section 3, we first recall 
the scaling behavior observed when the noise is colored and show explicitly that the standard averaging method 
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fails in this case. We then develop an averaging method that allows us to deal with an Ornstein-Uhlenbeck noise. 
This technique is based on a recursive transformation of coordinates where the noise itself is treated as a dynamical 
variable. An effective Langevin dynamics is derived for the slow variable that leads to an analytic expression for the 
P.D.F. and allows us to calculate the behavior of physical observables in the long time limit. Our results are finally 
extended to include the effect of a small dissipation in the system. 

II. REVIEW OF EARLIER RESULTS 

In this section, we review previously derived results that will be relevant for our present work. 

A. The model 

We consider a nonlinear oscillator of amplitude x(t), trapped in a confining potential U(x) and subject to a 
multiplicative noise £(£) : 

£«<)~^U «<)£«). (1) 

The statistical properties of the random function £(t) need not be specified at this stage. We restrict our analysis to 
the case where the potential hi behaves as an even polynomial in x when |ar| — ■> oo. A suitable rescaling of x allows us 
to write 

x 2n 

U with n > 2 . (2) 

2n 

As the amplitude x(t) of the oscillator grows with time, the behavior of U{x) for |a;| — > oo only is relevant and Eq.QJ 
reduces to 

^x{t)+x{tf n -^x{t)m- (3) 

B. Energy-angle coordinates 

An important feature of Eq.@ is that the deterministic system underlying it (obtained by setting £ = 0) is 
integrable. The associated energy and the angle variable are given by [Tll | 



, /E l/2n 



du 



E = + ~n~ n > and - / , , (4) 
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where the angle <f> is defined modulo the oscillation period 4:K n , with 
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K n = V^ n (5) 
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In terms of the energy-angle coordinates (E, (f>), the original variables (x, x) are given by 

x = E^ 2n S n (<!>), (6) 
x = {2n)^E 1 ' 2 SM , (7) 
where the hyperelliptic function S n is defined as 

V 2n 



The function S n and its derivative with respect to 4>, S' n , satisfy the following relation 



The presence of external noise spoils the integrability of the dynamical system © but does not preclude the use of 
(E, </)) instead of (x, i) as coordinates in phase space. Introducing an auxiliary variable f2 defined as 

Q = (2n) s ^ L , (10) 

equation J3J is written as a system of two coupled stochastic differential equations [ill Il2j | 

n = (n-l)S n (cl>)S' n (4>)m, (11) 
; n S n (cpf 
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This system is rigorously equivalent to the original problem (Eq. and has been derived without any hypothesis on 
the nature of the parametric perturbation which can even be a deterministic function or may assume arbitrary 
statistical properties. This external perturbation £(t) continuously injects energy into the system. We have shown 
analytically in [Til ll2T| that, when £ is a Gaussian white noise or a dichotomous Poisson noise, the typical value of Q 
grows algebraically with time. We also verified numerically that the same behavior is true for an Ornstein-Uhlenbeck 
noise. Therefore, as seen from Eq. I|12f> . the phase (f> is a fast variable in all cases of interest. Assuming that, in the long 
time limit, <f) is uniformly distributed over the interval [0, AK n ] of a period, we obtain, after averaging Eqs. © and l[7jl 
over the angle variable, the following equipartition relations 

<£} = ^<i 2 >, (13) 
2n 

(i 2 ) = (x 2n ). (14) 
These identities are in agreement with numerical simulations [llj |. 

C. Asymptotic formula for the P.D.F. in the white noise case 

We recall here the results obtained in [ll| for the case where £ (t) is a Gaussian white noise of zero mean value and 
amplitude D: 

m = o, 

mW)) = V5(t-t>). (15) 

Integrating out the angular variable <f> from the Fokker-Planck equation for the P.D.F. P t (Cl,(p) associated with the 
system <|11I12|1 . an averaged Fokker-Planck equation for the marginal distribution Pt(Q) is derived ^lj 



r(.) being the Euler Gamma function |l4 

We define, for sake of conciseness, the following two parameters 
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^ = Sl = (2n)n ) > ) ' , (17) 
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— T , . 2 M 2n / V 2n , 
M4 = S* = (2n)n ) - ( ) r ( , (18) 
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where the overline denotes the average over an angular period, for example S^ L — j^- Jq K " 5^(0)d</> . The relations Ijl7|) 
and H18(l are obtained by making the change of variable u = S n (4>), using Eq. © and evaluating the Eulerian integral 
of the first kind (i.e., Beta function) thus obtained [TlL[T^ |. 



Using Eqs. 1(16(1 and 1(17(1 . we deduce that the Fokker-Planck equation (|16(l corresponds to the following effective 
Langevin dynamics for the slow variable CI 

• , , i (n — 1)2? «2 (2n)^ (n — 1) , , . , . 

fi = (2n s i— i- % + V ; J_ >- VJ^V(t) , 19 
n + 3 11 y/n ,+ 3 

where 77(f) is an effective Gaussian white noise with amplitude T>. 

1 

From Eq. (|19|l . it is clear that the variable £1 has normal diffusive behavior with time: CI ~ (2?f) 2 . The averaged 
Fokker-Planck equation 1(16(1 is exactly solvable, leading to the following expression for the energy P.D.F. 
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Thus, the asymptotic time dependence of all moments of the energy, amplitude or velocity can be calculated analyt- 
ically and agree with numerical simulations [llj. In particular, using Eqs. I)l()fl. 113(1 and ((14(1 the following scaling 
relations are obtained 

E ~ (Dt)^, 

x - (Dt)^^ . (21) 

The physical observables grow algebraically with time, and the associated anomalous diffusion exponents depend only 
on the behavior of the confining potential at infinity. 

III. AVERAGING METHOD FOR COLORED NOISE AND CALCULATION OF THE P.D.F. 

We now consider £(f) to be a colored Gaussian noise with correlation time r that is obtained from the Ornstein- 
Uhlcnbeck equation 

^ = ->->. 

where 77(f) is a Gaussian white noise of zero mean value and of amplitude T>. In the stationary limit, when t, t' 3> t, 
we find: 

(£(*)) = and mW)) = ^ e-l*"*'!/- . (23) 

In section lill Al we recall qualitative scaling results derived in |12| for colored noise and show that a straightforward 
elimination of the angular variable 0, along the lines described above for white noise, leads to erroneous results. 
We then elaborate an averaging scheme that allows us to derive the long time behavior of the nonlinear oscillator 
subject to multiplicative Ornstcin-Uhlcnbcck noise : in subsection IHI Bl we define a new set of variables on the 
three dimensional phase space {Cl,4>, £) that retains the relevant correlations between the noise and the fast angular 
variable. In subsections IIII Cl and IIII Dl we derive the averaged Fokker-Planck equation and the associated effective 
low-dimensional Langevin system, respectively. Analytical results are obtained in IIII El for the non-dissipative system 
and extended in IIII Fl for a small dissipation rate. 

A. Qualitative behavior in the presence of colored noise 

In |l2j| . we deduced from a self-consistent scaling Ansatz that, when the noise has a finite correlation time r, the 
variable Cl has a subdiffusive behavior with time and scales as 



ci - (vty 



(24) 



Using this equation and Eqs. IjlOfl. (|13[1 and (|14|) . the following scaling relations are obtained 
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(25) 



where the factor Vt/2r 2 is found by dimensional analysis. Thus, the anomalous diffusion exponents are halved when 
the noise is colored. These colored noise scalings are observed when t — ► oo, even if the correlation time r is arbitrarily 
small. More precisely, the crossover between the white noise scalings (12111 and the colored noise scalings (I25|l occurs 
when the period T of the underlying deterministic oscillator (which is a decreasing function of its amplitude) is of the 
order of r, i.e., for a typical time t c ~ (Vt 2 )^ 1 . When t -C t c , the angular period of the system is much larger than 
the correlation time of the noise, which thus acts as if it were white. When t 3> t c , the noise is highly correlated over 
a period and its effect is smeared out leading to a slower diffusion. 

The qualitative scalings (|25(l have been obtained by elementary arguments [TTL Il2j|. However, an analytical calcu- 
lation of the P.D.F., that would yield quantitative formulae for the physical observables in the long time limit, has 
remained out of reach. Indeed, the standard averaging technique, which was successfully applied to white noise, fails 
for the Ornstein-Uhlenbeck process, as we show in the following. 

The random oscillator (fllll 2|) and the Ornstein-Uhlenbeck process Q22JI form a three dimensional stochastic system 
driven by a white noise r)(t). The Fokker-Planck equation for the joint P.D.F. Pt(0, is given by 

We now average this Fokker-Planck equation over the angular variable <p assuming that the probability measure for 
(f) is uniform over the interval [0, 4K n ] when t — > oo. We use the fact that the average of the derivative with respect 
to <f) of any function is zero: 



^(...) = 0. (27) 

This implies in particular that 



S n {<t>)S' n {<} ) ) = -d^Sm = 0. (28) 



Using these properties, we obtain the evolution equation for the marginal distribution P t (f2,£) 

dP t ldZPt V d 2 P t 



(29) 



dt t d£ 2t 2 d£ 2 ' 

This phase-averaged Fokker-Planck equation corresponds to the following Langevin dynamics for the variables (O, £) 

O = 

£ = -k--v(t)- (30) 

T T 

This result predicts that f2 is not stochastic anymore and is conserved. The integration over the angular variable 
averages out the noise itself and leads to conclusions that are blatantly wrong. The reason for this failure is the 
following: when the period of the underlying deterministic oscillator is less than r, the angular variable <p becomes 
fast as compared to both the energy E and the noise £. The straightforward elimination of <fi disregards the correlation 
between cf> and £ and eliminates the noise as well. A correct averaging scheme that takes into account the correlation 
between cj> and the noise £ will be developed below. 



B. Transformation of the equations 



We shall define recursively a new set of variables on the (f2, <f>, £) phase space that permits the adiabatic elimination 
of the angular variable without eliminating the relevant correlations between (f> and the noise £. Defining 



we deduce the dynamics of Y from Eqs. and i(12() 

r = (n - 1) — 5raW 2 € 4 = 2(n - l)S n {<p)S n {<f))^e, + + 0(Y ' ) , (32) 

1 F 

we have neglected here terms smaller than Y^ 1 (i.e., of order strictly higher than Y^ 1 ). The terms S n (4>)S' n (4>)(f> and 
S^(cj))S' n ((j))(f) being exact derivatives with respect to time, we integrate Eq. l|52"j l by parts and use Eq. (|2"2|) to obtain 

^-y = + °(y- 3/2 ) (33) 

Collecting all time derivatives on the left hand side (l.h.s.), we rewrite Eq. 1(33() as 

^ - S^K - = w| + (<») - /*) f + J + sto)f£±* + o(r- 3/2 ) , (34) 

where /12 is defined in Eq. (|17|l . The function (<S 2 (0) — ^2) is periodic in and its angular average vanishes identically 
by virtue of Eq. 1)17(1 . However, the product of this term with the noise £, that appears in Eq. 1(34(1 . does not average 
to zero (because of correlations between the angle variable and the noise). Therefore, defining a 0-periodic function 
T> n ((j)) that satisfies 

dV n {4>) 



=Si((f)-^ and 2? n (0)=O, (35) 
we rewrite the product (S^((f)) — ^2)^ as follows 

, c 2,^ u d^ P„(flg \ X> rl (0)(g + 77) S^)- /i 2 ^(0) + (n- l)V n {4>)S n {4>)S' n (0) ^ 2 nrv - 3 /2^ 

This relation can be verified by calculating the time derivative that appears on the right hand side (r.h.s.) and 
neglecting all terms of the order 0(y -3 / 2 ). Substituting Eq. 1(36(1 in Eq. 1(34(1 we finally obtain 



dT l^T r - s - m - T(2 n)-^yi/2 " 



$ PnWe 2S*(0) - ^(0) + (n - l)P n (0)S n (0)SM) 2 

M2 r r 2 (2n)-^yV2 T y 5 

^ r(2n) 2„y!/2 Y )t 

Considering the l.h.s. of this equation, it is natural to define a new dynamical variable Z such that 

Z = Y — (n — Dfam + ^ - + . (38) 

The two variables y and Z are identical at leading order. The dynamical equation for Z is obtained by writing Y 
in terms of of Z in Eq. 1(37(1 and neglecting all terms of order strictly higher than Z . This change of variable is 
straightforward on the l.h.s. of Eq. 1(37(1 and, because Y appears only in the denominator on the r.h.s., it is legitimate 
to replace Y by Z (at order Z^ 1 ). Finally, we obtain the following system of coupled Langevin equations 

1 -Z = J z (Z,ct>,i)+V z {Z,ct>,0—, (39) 



n — 1 r 

= J<f,(Z,cl>,0, (40) 

£ = --C--ry(i), (41) 
r r 



where the current and diffusion functions are given by 

J Z {Z,4>,0 = M2 - i i ^W-^W + t"- 1 ^^-^^) ^ (42) 

^>z{Z,^i) = S n {(f>) + _L^ 1/9 + 7 > ( 43 ) 

T{2n) 2 n 'Z i -i z ^ 

J<l>{Z,(p,Z) = — —j £. (44) 

(2n)» S2(Z,0,£) 

We now write the Fokker-Planck equation associated with this stochastic system and average it over the rapid 
variations of the angular variable <j). 



C. Averaging the Fokker-Planck equation 

The Fokker-Planck equation for the P.D.F. H t (Z,<j),£) corresponding to the system ((391401 and I41J1 is given by 

~dt = - {n - l) dz {JzI{t) -d4> {J ^ t)+ T^r + 

- v£f>'m*>™ ~ {n ~ ~ {n - l) + P} ■ (45) 

We now integrate out the fast angular variable (j> from this equation, i.e., we consider that in the long time limit, the 
P.D.F. Ht(Z, <f>, £) becomes uniform in <fi and reduces to the function flt(Z, £) of the two variables Z, £. Because the 
stochastic system 1(391 1401 and I41|) is valid up to the order we should retain, while averaging Eq. I(45|l. only 

the terms that scale at most as Z -2 (we recall that the derivative operator itself scales as Z^ 1 ). 

We now calculate the angular average of each term appearing on the r.h.s. of Eq. 145(1 . To deal with the first term, 
we need to calculate the average of the current function defined in Eq. 1(42(1 and we obtain 



1-(7 e\ VJM , 254(0) - m Sl{$) + (n- l)P n (0)S n (0R(0) rf 

JzK/'yZ,) = M2- + — - — r— tt H ~ ? ■ (4b) 

t T 2 (2n)"2^Z 1 /2 tZ 

Integrating by parts and using Eq. 1(35(1 . we write 

Vj^mM) = -^. 2 (0)(W) - M2 ) = ^i^i , (47) 
where /i4 was defined in Eq. {TSJ. Using Eqs. JTZJ, {ISJl I^SJ) and l|I7|l. Jz(Z,£) is written as 

^,0=M2;++^(2/i4-^ + ^(M2-M4)) ■ (48) 
The average of the diffusion function defined in Eq. ((43(1 can also be calculated in a similar manner : 

P^,£)= M2 + M. (49 ) 
From this equation the angular averages of the second derivative terms that appear in Eq. 1(45(1 are readily found 



jl^Adw) = ,50) 



az 2 ^ = az( (Ai2 + — } ^J- (52) 

In Eq. 1(50(1 we have retained only the leading order in the average of the square of the diffusion function T>z, because 
the second derivative with respect to Z already scales as Z~ 2 . 



Finally, we deduce from Eq. (|27|l that 



^(J^n t )=0. (53) 



Inserting the expressions obtained in Eqs. 153|) . (|48|) . J5UJ), i(51|) and 1(521) in Eq. I|45|) . we complete the derivation of 
the averaged Fokker-Planck equation : 

-5r = -( n -!)^ — £ + ^ C n 



^ v J dZ\ \ r 2rZ s / t d£ 

2^|(« - D Vg^ - (» - 1)^((M2 + f f)n«) - (» - D^((M 2 + + ^ j ■ (54) 

This equation describes a stochastic motion in the two-dimensional phase space (Z,£). We shall now write the effective 
Langevin equations for (Z, £) that correspond to this averaged Fokker-Planck equation. 



D. Effective Langevin equations 



Consider the following stochastic system 



n — 1 t 2tZ r V Z / r 

£ = --e --»»(*). (56) 
r r 

where ?7i(t) and r)2(t) are two independent Gaussian white noises of amplitude T>. The Fokker-Planck equation 
associated with the stochastic system 1)551 156|) is identical to Eq. (|54|l but for the second derivative in Z. In Eq. I|54|) 
this term is given by 

whereas the corresponding term for the Fokker-Planck equation corresponding to Eqs. (|55l I56J) is 

JL (n _ 1} J (w _ + £ („ + 1 { ) ^ (, 2 + ^ n,} . , 58) 



Retaining in Eq. I|58|) only contributions of order up to Z~ 2 , as done earlier, we observe that the two expressions 
coincide. Hence, the Fokker-Planck equation corresponding to the system (|55I56|) is precisely given by the averaged 
Fokker-Planck equation l|54(l in the two-dimensional phase space (Z, £) obtained after adiabatic elimination of the fast 
angular variable. This Langevin dynamics for the slow variable Z plays a role similar to that of Eq. 1)19(1 in the case 
of white noise. However, the stochastic system ((55156(1 does not admit an exact solution and we need to make further 
simplifications in order to obtain analytic results. The main idea is to eliminate the white noise 772 (i) from Eq. 155(1 
so that the equation for Z is partially decoupled from that for £. We therefore define a new variable Z2 such that 



Z 2 = Z+(n-l)^ + (n-l)^e- (59) 



The time evolution of Z2 is given by 



^ = (n-l)(3-n)(M4-/ij) ^ y/^-^ 
2tZ 2 t 

To derive this equation we retained only the terms that scale at most as Z^ 1 ■ Using the fact that the mean value 
(£ 2 ) of the Ornstein-Uhlenbeck process is equal to T>/2t, we write Eq. ((60(1 as follows 



(n-l)(3-n)£>M4-M2 , , , n V/M^mI , a , (n~ 1)(3 - n){^ - Mp /xa /t 2\\ /fin 
4t z Z 2 t 2rZn 



This equation contains two random noise sources: the white noise rji(t), and the nonlinear colored noise £ 2 — (£ 2 ), of 
zero mean and of finite variance. This colored noise is multiplied by a prefactor I/Z2 and, therefore, its contribution 
in the long time limit becomes negligible as compared to that of the white noise. We thus discard this term and 
obtain the following white noise effective Langevin dynamics for Z2 



^ = (n - 1)(3 - n)V fi 4 - f4 + t n _ x n - l4 . 



4r2 Z2 An-D^^mit). (62) 

This equation has the same mathematical structure as Eq. (|19H obtained for multiplicative white noise, but the 
coefficients in these equations are different. Besides, although the effective variables and Z2 used, respectively, 
for white noise and colored noise, satisfy similar equations, they are not equivalent [in fact, from Eqs. 1)381 and 13 1|) 
we observe that Z scales as J7 1 / 2 ]. Thus, the difference between the white noise and the colored noise problems is 
embodied in the successive transformations that relate Q and Z2 to the original energy-angle coordinates. 



E. Analytical results 

We now derive the analytic expression for the P.D.F. of the energy of the system in the long time limit and 
then calculate the asymptotic behavior of physical observables such as mean energy, mean position-square and mean 
velocity square. 

Observing that the Fokker-Planck equation associated with the effective Langevin dynamics of Z2, given by Eq. I|62p. 
is exactly solvable, we deduce the following expression for the P.D.F. 

Pti Z 2 ) = , 2 . J- f#) ^ exp f - , (63) 

v ' r f n+i\ ,/2Af,\2AtJ l \ 2 At/ ' v ; 



( n+1 \ 



with 



A = (n- 1) 2 (^4 -A»a)^- ( 64 ) 

The expression for the energy E as a function of Z2 can be derived from Eqs. O, |(3TJ, J2U and g2J. At leading 
order, we obtain 

Z 2 = 2nE^ , (65) 

this relation is valid when £ > 1, i.e., in the long time limit. We thus deduce the asymptotic expression, valid for 
t — ► 00, of the probability distribution function of the energy : 

Pt(E) = nT ( J\ (W ^ 6XP ( At ) ' ^ 

If we compare this expression of the P.D.F. for colored noise with the one derived for white noise H20|) . we observe 
that the two distribution functions are of the type E a exp (—cE^/Dtj, with the characteristic exponents a and f3 and 
the constant c depending precisely on nature of the noise. 

The P.D.F. (|66() together with the assumption of uniform angular measure in the asymptotic regime leads to 
analytical expressions for the mean value of any observable of the system. In particular, from Eqs. ©, l|13|l and 166() 
we derive the statistical mean of the position, velocity and energy of the system in the long time limit 

ITA - F (jj^jy) / At \ _ F (jjjfg) / (n - 1) 2 (A^4 - t4)Vt \ 2(^T) . 

K ) v ( n+i \\2n 2 ) T ( n+i \ V 2n 2 r 2 ) ' 1 ' 

1 \M.n-l)J \*(n-l)J 



a 2 ) = 



'In 



(E) , (68) 
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FIG. 1: Mean value of: (a) the energy and (b) the square amplitude of the non-dissipative nonlinear oscillator subject 
to parametric Ornstein-Uhlenbeck noise, with T> = 1 and r = 5. The insets show the following ratios vs. time t: (a): 
(E)/{Vt/T 2 ) n,( ' Kn ~ 1)) \ (b): (x 2 >/(©t/r 2 ) 1/(2( "~ 1)) . The horizontal dotted lines correspond to the numerical prefactors given 

in Eqs. mmm\ . 



The power-law behaviors predicted by these analytical expressions are in accordance with the qualitative scaling 
relations of Eq. deduced from heuristic arguments in [l^]. Giving some particular values to the parameter n, we 



obtain 



for n = 2, 



, (E) = 0.0467 (^) , (i 2 ) = 0.0623 (^) , (x 2 ) = 0.169 (^) 



<T>t\V* 
7z 



(70) 



'T> t\ 3 / 4 

for n = 3, (E) = 0.0835 11 j 



, (x') = 0.125 ( —J , {a?) = 0.296 ^ , 



(71) 



for n = 4, (£) = 0.0933 — 



Vt\V3 



= 0.149 



(x z ) = 0.338 



2?t\l/6 



IT 5 " 



This results are in excellent agreement with numerical simulations (see Figure 

F. Extension to the case with small dissipation 



(72) 



The quantitative analysis described above can be generalized to the study of a stochastic nonlinear oscillator with 
very small dissipation. Introducing a linear friction term with dissipation rate 7 in Eq. we obtain 



±. x{ t )+1 ± x{ t)+ X (t)^ = x(t)m- 



(73) 



where £ is an Ornstein-Uhlenbeck noise. This second order random differential equation can be written as a stochastic 
system in the (O, 4>) coordinates |l2j 



n = _ 7 _- Tl s;(0) 2 fi + (n-i)5„(0)5;(^^), 

(2n)» 

• s n (<t>)s> n (<t>) , n 



(2n)» (2n)" 



(74) 
(75) 



Saturation is attained when the power injected by the random force £(i) is balanced by dissipation. The variable 
then reaches a finite mean-value and fluctuates around this value. The separation between fast and slow variables 



is valid only in the limit of a vanishingly small dissipation rate : in this case, fl saturates to a very large value, the 
term 0/(2n)« is the dominant term in Eq. (|75|l . and therefore the phase <f> varies rapidly with time. We thus derive 
an analytic expression for the stationary P.D.F. of the stochastic process (|73|) under the hypothesis 7 <C 1, i.e., in 
dimensionless units, we assume that 72?~ 1 / 3 <C 1 and 71 < 1. It is then possible to repeat all the calculations of 
sections UTTbI InTcl and IHTDl in the presence of a dissipative contribution and to derive the following effective white 
noise Langevin dynamics for the variable Z2, defined in Eq. (|59|l 



7 o n ~ 1 7 (n - - n)V m - $ VM4 - Ml , a , 7R , 

2 7 ^Ti 2 ( } — t — m ®' ( ' 

The stationary Fokker-Planck equation corresponding to this dynamics can be explicitly solved and the following 
expression for the stationary P.D.F. of the energy is obtained 

. . 2(n-l) /2n 2 \4r^i7 i-n / 2n 2 E 2 ^\ , , 

P ^ {E) = nr ( J) M ^ exp ( S— )' (77) 



where A is given by 



A -s^ny A -<-"- 1 >(«-'«j^- (78) 



We emphasize that the P.D.F. given in Eq. (|77|l is not of the canonical Gibbs-Boltzmann form and therefore does 
not represent a state of thermodynamic equilibrium. Using this P.D.F., the moments of the energy, position and 
velocity can be calculated. The analytical expressions of (E), (x 2 ) and (x 2 ) can be deduced in a formal way from 
Eqs. H67JI. <|68[1 and l|69|) by substituting for the time variable t the expression 

1 n + 1 ( 79 ) 
7 4(n-l) 1 1 

This expression defines in fact a dissipative time scale td such that for t <C td the system behaves as if it were non- 
dissipative and for t 3> td, the system is settled in a non-equilibrium stationary state. For t ~ td, the time-dependent 
P.D.F. I|66|) matches the stationary P.D.F. 1|77[1. We emphasize that an important assumption in our calculations is 
that the P.D.F. is uniform with respect to the fast angular variable in the long time limit. This hypothesis breaks 
down when the damping rate 7 is high such that the angle and the energy vary on comparable time scales. The 
separation between fast and slow variables is then no more possible. In fact, when 7 exceeds a critical value (keeping 
all other parameters fixed) the system undergoes a noise-induced phase transition and the stationary state reduces to 
the fixed point x = x = 0. Although for white noise this critical value is exactly known, an exact calculation of 
this bifurcation threshold for colored noise remains open question [16) , The recursive adiabatic elimination technique 
developed here may perhaps be useful to tackle this challenging problem. 



IV. CONCLUSION 



We have derived in this work analytical results for the nonlinear oscillator subject to a multiplicative colored noise. 
In the problem we have studied here, the following features play an essential role : 

(i) The model's equation is genuinely of the second order because we have considered the case of zero (or very 
small) damping. 

(ii) The dynamics is necessarily nonlinear because we have studied the behavior at large times, for which the 
nonlinearity of the confining potential becomes relevant. 

(iii) In the long-time limit, the finite correlation time r of the multiplicative Ornstein-Uhlcnbeck process is not the 
smallest time scale of the problem (the angular period defines a smaller time scale) . Therefore, the noise is genuinely 
colored and can not be treated as a 'quasi-white' even if t is arbitrarily small. 

We emphasize that analytical results for such second-order Langevin equations are scarce. Moreover, in previous 
works the above features are not present simultaneously. Indeed, different approximations (for a review see, e.g. , 
[TtI ^|) necessarily neglect one or more aspects of the problem. For example, in the large damping limit [T^ |. 
the inertial term is neglected and the model is reduced to a first order Langevin equation. Short-time behavior is 
dominated by the quadratic part of the confining potential, thus leading to a linear stochastic equation. In small 
correlation-time expansions that provide effective Fokker-Planck equations [20I l2l| , it is implicitly assumed that r is 



the smallest time scale; but perturbative expansions in the vicinity of the white noise limit are unable to predict the 
anomalous diffusion exponents for colored noise. Specific techniques such as the Unified Colored Noise Approximation 
|22l l2.'"'J| appropriate for additive noise problems, seem to be unsuitable for multiplicative noise with small dissipation. 

Thanks to a recursive adiabatic averaging method, we have been able to perform a mathematical analysis of the 
long time behavior of a nonlinear oscillator with parametric noise in a non-perturbative manner. We have obtained 
an explicit formula for the P.D.F. in the long time limit which allows us to calculate the statistical average of any 
physical observable. Our results have been successfully compared to numerical simulations; the agreement is not only 
qualitative but also quantitative. In the case of a dissipative system, we have derived an expression for the stationary 
P.D.F. that compares well with numerical simulations in the weak damping regime. 

This study concludes a series of works in which we have generalized the adiabatic averaging technique 

to derive analytical results for quasi-Hamiltonian nonlinear random oscillators subject to an additive/multiplicative, 
white or colored Gaussian noise. 
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